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RANDOM POLYNOMIALS OVER R AND OVER C. 



DIEGO ARMENTANO 

Abstract. This article is divided in two parts. In the first part 
we review some recent results concerning the expected number of 
real roots of random system of polynomial equations. In the sec- 
ond part we deal with a different problem, namely, the distribution 
of the roots of certain complex random polynomials. We discuss 
a recent result in this direction, which shows that the associated 
points in the sphere (via the stereographic projection) are surpris- 
ingly well-suited with respect to the minimal logarithmic energy 
on the sphere. 



1. Introduction 

Let us consider a system of m polynomial equations in m unknowns 
over a field K, 

(1) f t (x):= J2 a f xj (i = l,---,m). 

Ib'll<* 

The notation in (TjQ) is the following: x := (x±, . . . , x m ) denotes a point 
in K m , j := . . . ,j m ) a multi-index of non-negative integers, ||j|| = 

Y^h=i3hi x ^ — x ^ ' ' -& m , Cij = aj'' r .j m , and d{ is the degree of the 
polynomial fa. 

We are interested in the solutions of the system of equations 

(2) fi(x) = (i = l,...,m), 

lying in some subset V of K m . Throughout this review we are mainly 
concerned with the case IK = R or IK = C. 

If we choose at random the coefficients {a^}, then the solution of the 
system (j2J) becomes a random subset of lK m . This is the main object 
of this review. 

In the first part of this paper we focus on the real case. The main 
problem we consider is that of understanding N*(V): the number of 
solutions lying in the Borel subset V of lR m . 
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In the second part we deal with a different problem: How are the 
roots of complex polynomials distributed? 

This article is organized as follows: 
In Section [H we start with some historical remarks on random polyno- 
mials. After that we move to the case of random systems of equations. 
We mention some recent results for centered Gaussian distributions. 
In Section \2J\ we consider the non-centered case, which has also been 
called "smooth-analysis" in the last years. That is, we start with a fixed 
(non-random) polynomial system, then we perturb it with a polyno- 
mial noise, and we ask what can be said about the number of roots of 
the perturbed system. In Section we review a result which com- 
putes the expected number of roots of a random system of polynomial 
equations expressed in a different basis, namely, the Bernstein basis. 
Finally in Section&we focus on the complex case. We discuss a recent 
result concerning the distribution of points in the sphere associated 
with roots of random complex polynomials. 

This review follows the talk given by the author in the colloquium 
which was held the inauguration of the Franco- Uruguayan Institute of 
Mathematics, in Punta del Este, Uruguay, on December 2009. 

2. The Number of Real Roots of Random Polynomials 

The study of the expectation of the number of real roots of a random 
polynomial started in the thirties with the work of Block and Polya 
[7]. Further investigations were made by Littlewood and Offord [13] . 
However, the first sharp result is due to M. Kac (see Kac [TTl fT2"] ). who 
gives the asymptotic value 



when the coefficients of the degree d univariate polynomial / are Gauss- 
ian centered independent random variables N(0, 1) (see the book by 
Bharucha-Reid and Sambandham [6]). 

The first important result in the study of real roots of random system 
of polynomial equations is due to Shub and Smale [2U] in 1992, where 
the authors computed the expectation of A r ^(M m ) when the coefficients 
are Gaussian centered independent random variables having variances: 



E(N f (R)) « -logd, as d ->■ +oo, 



(3) 




dil 



Their result was 



(4) 



E(N f (R m )) = yfd x ■ ■ ■ d 
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that is, the square root of the Bezout number associated to the sys- 
tem. The proof is based on a double fibration manipulation of the 
co-area formula. Some extensions of their work, including new re- 
sults for one polynomial in one variable, can be found in Edelman- 
Kostlan[10]. There are also other extensions to multi-homogeneous 
systems in McLennan [16], and, partially, to sparse systems in Rojas[T7] 
and Malajovich-Rojas|15j. A similar question for the number of critical 
points of real- valued polynomial random functions has been considered 
in Dedieu-MalajovichJH]. 

The probability law of the Shub-Smale model defined in §5$) has 
the simplifying property of being invariant under the action of the 
orthogonal group in M m . In Kostlan[13] one can find the classification 
of all Gaussian probability distributions over the coefficients with this 
geometric invariant property. 

In 2005, Azai's and Wschebor gave a new and deep insight to this 
problem. The key point is using the Rice formula for random Gauss- 
ian fields (cf. Azai's- Wscheb or jS]). This formula allows one to extend 
the Shub-Smale result to other probability distributions over the coef- 
ficients. A general formula for K(N*(V)) when the random functions 
fi (i = 1, . . . ,m) are stochastically independent and their law is cen- 
tered and invariant under the orthogonal group on M m can be found 
in Azai's- Wscheb or [1]. This includes the Shub-Smale formula @ as 
a special case. Moreover, Rice formula appears to be the instrument 
to consider a major problem in the subject which is to find the as- 
ymptotic distribution of AT* (V) (under some normalization). The only 
published results of which the author is aware concern asymptotic vari- 
ances as m — > +oo. (See Wschebor[2S] for a detailed description in this 
direction and a simpler proof of Shub-Smale result). 

2.1. Non-centered Systems. The aim of this section is to remove 
the hypothesis that the coefficients have zero expectation. 

One way to look at this problem is to start with a non-random system 
of equations (the "signal" ) 

(5) P*(x) = (i = l,...,m), 

perturb it with a polynomial noise X^x) (i = 1, . . . , m), that is, con- 
sider 

Pi(x) + Xi{x) = (z = l,...,m), 

and ask what one can say about the number of roots of the new system, 
or, how much the noise modifies the number of roots of the determin- 
istic part. (For short, we denote N* = N^(M. m )). 
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Roughly speaking, we prove in Theorem U\ that if the relation signal 
over noise is neither too big nor too small, in a sense that will be made 
precise later on, there exist positive constants C, 9, where < 9 < 1, 
such that 



(6) 



< C9 m E(N 



x \ 



Inequality (|6]) becomes of interest if the starting non-random system 
(JHJ) has a large number of roots, possibly infinite, and m is large. In 
this situation, the effect of adding polynomial noise is a reduction at 
a geometric rate of the expected number of roots, as compared to the 
centered case in which all the Pj's are identically zero. 

For simplicity we assume that the polynomial noise X has the Shub- 
Smale distribution. However, one should keep in mind that the re- 
sult can be extended to other orthogonally invariant distributions (cf. 
Ar ment ano-Wscheb or [2] ) . 

Before the statement of TheoremU\below , we need to introduce some 
additional notations. 

In this simplified situation, one only needs hypotheses concerning 
the relation between the signal P and the Shub-Smale noise X, which 
roughly speaking should neither be too small nor too big. 

Since X has the Shub-Smale distribution, from fl3]) we get 



Var(X,(x)) 
Define 



1 



2\di 



K(P t ) := 



sup 



Vx e 



m). 



P 



sup 

xeM m \{0} 



■1' 



1 + ||x|| 2 )^/ 2 
d_ ( P 
dp 



for i 



m, where 



is the Euclidean norm, and S- denotes the 



derivative in the direction defined by 
For r > 0, put: 

P l (x) 2 



L(Pi,r) 



inf 



\2\di 



dp 

at each point i^O. 



(i = 1, . . . ,m). 



IWI>r (1 + \\X\ > 

One can check by means of elementary computations that for each P 
as above, one has 

H{P) < oo, K(P) < oo. 

With these notations, we introduce the following hypotheses on the 
systems as m grows: 
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Hi) 



(7a) 



1 m 

= -• V 

i=l 




o(l) as m — >■ +00 



(7b) 



B, 



'in 




K 2 



1 



(Pi) 



o(l) as m — > +00. 



H2) There exist positive constants tq, £ such that if r > tq: 
L(Pi, r) > I for all % — 1, . . . , m. 
Theorem 1. Under the hypotheses Hi) and H2), one has 



where C, 6 are positive constants, < 6 < 1. 

2.1.1. Remarks on the statement of Theorem^ 

• It is obvious that our problem does not depend on the order in 
which the equations 



appear. However, conditions fl7a|) and (17b"]) in hypothesis H 3 ) 
do depend on the order. One can state them by saying that 
there exists an order % — 1, . . . , m on the equations, such that 
(17a]) and flSjJ) hold true. 

• Condition Hi) can be interpreted as a bound on the quotient 
signal over noise. In fact, it concerns the gradient of this quo- 
tient. In ( 17b |) the radial derivative appears, which happens to 
decrease faster as — > 00 than the other components of the 
gradient. 

Clearly, if H(Pi), K(Pj) are bounded by fixed constants, (17a]) 
and ( 17b |) are verified. Also, some of them may grow as m — > +00 
provided ( 17aj) and ( |7b1) remain satisfied. 

• Hypothesis i/2) goes - in some sense - in the opposite direction: 
For large values of ||x|| we need a lower bound of the relation 
signal over noise. 

• A result of the type of Theorem [I] can not be obtained without 
putting some restrictions on the relation signal over noise. In 
fact, consider the system 



where a is a positive real parameter. If we let a +00, the 
relation signal over noise tends to zero and the expected number 
of roots will tend to E(N X ). On the other hand, if a I 0, E(N X ) 
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E(N P+X ) < C9 m E(N x ) 



Pi(x) +X l (x) =0 (i = l,...,m) 



(9) 



Pi(x) +aXi(x) = (i = l,...,m), 



(> 
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can have different behaviours. For example, if P is a "regular" 
system, the expected value of the number of roots of (JH]) tends 
to the number of roots of Pi(x) — 0, (i = 1, . . . , m), which may 
be much bigger than K(N X ). In this case, the relation signal 
over noise tends to infinity. 
• As it was mentioned before we can extend Theorem[J\to other 
orthogonally invariant distributions. However, for the general 
version we need to add more hypotheses. 

In the next paragraphs we are going to give two simple examples. 
For the proof of Theorem [7] and more examples with different noises 
see Armentano-Wschebor[2]. 

2.1.2. Some Examples. We assume that the degrees di are uniformly 
bounded. 

For the first example, let 



Pi(x) 



\x\\ di 



where di is even and r is positive and remains bounded as m varies. 
Then, one has: 



d_ 

dp 

V 



P 



(1+1^112)^/2 
Pi 



[X 



di \\x\\ d ' 1 + di r di \\x\ 

2^+1 



< 



di(l 



1 + \\x\\ 2 ) d >/ 2 
which implies 



(i + If.. 

di llxll^" 2 + di r d 



(1 + ||x|| 2 ) 3 / 2 



x 



Pi 



||a;||2)*/2 

Again, since the degrees di, . . 



[x 



< 



di{l + r di ) 



[1 + ||x|| 2 ) 3 / 2 ' 

d m are bounded by a constant that 
does not depend on m, Hi) follows. H 2 ) also holds under the same 
hypothesis. 



Notice that an interest in this choice of the Pj's lies in the fact that 
obviously the system Pi(x) — (i — 1, . . . ,m) has an infinite number 
of roots (all points in the sphere of radius r centered at the origin are 
solutions), but the expected number of roots of the perturbed system 
is geometrically smaller than the Shub-Smale expectation, when m is 
large. 

Our second example is the following: Let T be a polynomial of degree 
d in one variable that has d distinct real roots. Define: 



P i (x 1 ,...,x m ) = T(x i ) (i = l,...,m). 
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One can easily check that the system verifies our hypotheses, so that 
there exist C, 9 positive constants, < 9 < 1 such that 

E(N P+X ) < C9 m d m/2 , 

where we have used the Shub-Smale formula when the degrees are all 
the same. On the other hand, it is clear that N p = d m so that the 
diminishing effect of the noise on the number of roots can be observed. 
A number of variations of these examples for P can be constructed, 
but we will not pursue the subject here. 

2.2. Other Polynomial Basis. Up to now all probability measures 
were introduced in a particular basis, namely, the monomial basis 
{ x ' 7 }||ill<'i- However, in many situations, polynomial systems are ex- 
pressed in different basis, for example, orthogonal polynomials, har- 
monic polynomials, Bernstein polynomials, etc. So, it is a natural ques- 
tion to ask: What can be said about N* (V) when the randomization is 
performed in a different basis? 

For the case of random orthogonal polynomials see Barucha-Reid 
and Sambandham[6], and Edelman-Kostlan|10] for random harmonic 
polynomials. 

In this section following Armentano-Dedieu[5] we give an answer 
to the average number of real roots of a random system of equations 
expresed in the Bernstein basis. Let us be more precise: 

The Bernstein basis is given by: 

bdLkix) = f fV(l - *Y~ k i 0<k<d, 



k, 

in the case of univariate polynomials, and 

fe dJ (xi, . . . ,x m ) = (^j x i ■ ■ -z£r(l - X! - ... - x m ) d ~ M , \\j\\ < d, 

for polynomials in m variables, where j — (ji, ■ ■ ■ , j m ) is a multi- integer, 
and ( .) is the multinomial coefficient. 

Let us consider the set of real polynomial systems in m variables, 

fi(xi,...,x m )= cq'b d) j(xi, . . .,x m ) (i = 1, . . . ,m). 

\\j\\<di 

Take the coefficients to be independent Gaussian standard random 
variables. 
Define 

r : W n -»■ P (R m+1 ) 

by 

t(X\ , . . . , X m ) \X\ , . . . , X m , 1 X\ ... Im] • 
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Here P(K m+1 ) is the projective space associated with M m+1 , [y] is the 
class of the vector y G IR m+1 , y ^ 0, for the equivalence relation defining 
this projective space. The (unique) orthogonally invariant probability 
measure in P (M m+1 ) is denoted by A m . 

With the above notation the following theorem holds: 

Theorem 2. (1) For any Borel set V in M m we have 

E (N*{V)) = X m (r(V))y/d 1 ...d m . 

In particular 

(2) E (Nf) = Vdt.. d m , 

(3) E (A^(A m )) = ^Jd 1 ...d m /T\ where 

A m = {i6l ra : %i > and x x + . . . + x m < 1} , 

(4) When m = 1, for any interval I = [a, /3] C IR, one has 

E (N f (I)) = — (arctan(2/3 - 1) - arctan(2a - 1)) . 

The fourth assertion in Theorem [His deduced from the first assertion 
but it also can be derived from Crofton's formula (see for example 
Edelman-Kostlan [TU] ) . 

For the proof of Theorem{^see Armentano-Dedieu[3] 

3. Distribution of Complex Roots of Random Polynomials 

In this part we will see that points in the sphere associated with 
roots of Shub-Smale complex analogue random polynomials via the 
stereographic projection, are surprisingly well-suited with respect to 
the minimal logarithmic energy on the sphere. That is, they provide 
a fairly good approximation to a classical minimization problem over 
the sphere, namely, the Elliptic Fekete points problem. 

Next paragraphs follows closely Armentano-Beltran-ShubpQ, where 
one can find proofs and more detailed references. 

Given x±, . . . , x N G S 2 = {x e M 3 : \\x\\ = 1}, let 

(10) V(x 1 ,...,x N ) = hi Tl T 77 = - 5^ lnHxi-xJ 

l<i<j<N 11 1 l<i<j<N 

be the logarithmic energy of the iV-tuple x\, . . . , x^. Let 
V N = min V(x 1 ,...,x N ) 

denote the minimum of this function. iV-tuples minimizing the quan- 
tity (TIT)]) are usually called Elliptic Fekete Points. The problem of 



RANDOM POLYNOMIALS 



9 



finding (or even approximate) such optimal configurations is a classical 
problem (see White [23J for its origins). 

During the last decades this problem has attracted much attention, 
and the number of papers concerning it has grown amazingly. The 
reader may see Kuijlaars-Saff[19] for a nice survey. 

In the list of Smale's problems for the XXI Century [22], problem 
number 7 reads: 

Can one find Xi, . . . ,x^ G § 2 such that 

(11) V(x u ...,x N )-V N <c\nN, 

c a universal constant? 

More precisely, Smale demands a real number algorithm in the sense 
of Blum-Cucker-Shub-Smale[8j that with input N returns a iV-tuple 
X\, . . . , xpf satisfying equation pip , and such that the running time is 
polynomial on N. 

One of the main difficulties when dealing with this problem is that 
the value of Vn is not even known up to logarithmic precision. In 
Rakhmanov-S aff-Zhou [IB] the authors proved that if one defines Cn 
by 

V 2 , /4\ N\nN „ 

-0.112768770... < liminf C N < limsupCV < -0.0234973... 

Let X\, . . . , Xn be independent random variables with common uni- 
form distribution over the sphere. One can easily show that the ex- 
pected value of the function V(X 1 , . . . , X N ) in this case is, 

N 2 /4\ N 
(13) E(V(X h X N )) = -— In (- + — In 

Thus, this random choice of points in the sphere with independent 
uniform distribution already provides a reasonable approach to the 
minimal value Vn, accurate to the order of 0(N In N). 

On one side, this probability distribution has an important property, 
namely, invariance under the action of the orthogonal group on the 
sphere. However, on the other hand this probability distribution lacks 
on correlation between points. More precisely, in order to obtain well- 
suited configurations one needs some kind of repelling property between 
points, and in this direction independence is not favorable. Hence, 
it is a natural question whether other handy orthogonally invariant 



(12) 
then, 



V N = - 
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probability distributions may yield better expected values. Here is 
where complex random polynomials comes into account. 
Given z G C, let 
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g C x 



1 + 1*1 

be the associated points in the Riemann Sphere, i.e. the sphere of 
radius 1/2 centered at (0,0,1/2). Finally, let 

X = 25- (0,0,1) G § 2 

be the associated points in the unit sphere. 

Given a polynomial / in one complex variable of degree N, we con- 
sider the mapping 

f^V{X 1 ,...,X N ), 
where X^ (i — 1, . . . , N) are the associated roots of / in the unit sphere. 
Notice that this map is well defined in the sense that it does not depend 
on the way we choose to order the roots. 

Theorem 3. Let f(z) = ^2k=o akzk ^ e a complex random polynomial, 
such that the coefficients are independent complex random variables, 
such that the real and imaginary parts of a& are independent (real) 
Gaussian random variables centered at with variance (,). Then, 
with the notations above, 

v „ N\ /4\ N\nN N 4 
E(V(X ll ... l ^)) = -- r ln^J -_^ + T ln-. 

Comparing Theorem^ with equations (IT2"|) and (1131) . we see that the 
value of V is surpringsingly small at points coming from the solution 
set of this random polynomials. More precisely, necessarily many ran- 
dom realizations of the coefficients will produce values of V below the 
average and very close to Vn, possibly close enough to satisfy equation 

(ED- 

Notice that, taking the homogeneous counterpart of /, Theorem^ 
can be restated for random homogeneous polynomials and consider- 
ing its complex projective solutions, under the identification of IP(C 2 ) 
with the Riemann sphere. In this fashion, the induced probability 
distribution over the space of homogeneous polynomials in two com- 
plex variables corresponds to the classical unitarily invariant Hermitian 
structure of the respective space (see Blum-Cucker-Shub-Smale[8]). 
Therefore, the probability distribution of the roots in P(C 2 ) is invari- 
ant under the action of the unitary group. 

Is not difficult to prove that the unitary group action over P(C 2 ) 
correspond to the special orthogonal group of the unit sphere. Hence, 
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the distribution of the associated random roots on the sphere is orthog- 
onally invariant. Thus, Theorem [3] is another geometric confirmation 
of the repelling property of the roots of this Gaussian random polyno- 
mials. 

For a proof of Theorem and a more detailed discussion on this 
account see Armentano-Beltran-Shub[T]. See also Shub-Smale [21] . 
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